Chaotic, staggered and polarized dynamics in 
opinion forming: the contrarian effect 

Christian Borghesi 1 * and Serge Galam 2 ^ 
Service de Physique de PEtat Condens (CNRS URA 2464), 
CEA Saclay, 91191 Gif sur Yvette Cedex, France 
2 Centre de Recherche en Epistemologie Appliquee, 
CREA-Ecole Polytechnique (CNRS UMR 7656), 
1, rue Descartes, F-75005 Paris, France 

(15 Mar 2006) 



Abstract 

We revisit the no tie breaking 2-state Galam contrarian model of opin- 
ion dynamics for update groups of size 3. While the initial model assumes 
a constant density of contrarians a for both opinions, it now depends 
for each opinion on its global support. Proportionate contrarians are 
thus found to indeed preserve the former case main results. However, 
restricting the contrarian behavior to only the current collective major- 
ity, makes the dynamics more complex with novel features. For a density 
a < a c — 1/9 of one-sided contrarians, a chaotic basin is found in the 
fifty-fifty region separated from two majority-minority point attractors, 
one on each side. For 1/9 < a < 0.301 only the chaotic basin survives. In 
the range a > 0.301 the chaotic basin disappears and the majority starts 
to alternate between the two opinions with a staggered flow towards two 
point attractors. We then study the effect of both, decoupling the lo- 
cal update time sequence from the contrarian behavior activation, and a 
smoothing of the majority rule. A status quo driven bias for contrarian 
activation is also considered. Introduction of unsettled agents driven in 
the debate on a contrarian basis is shown to only shrink the chaotic basin. 
The model may shed light to recent apparent contradictory elections with 
on the one hand very tied results like in US in 2000 and in Germany in 
2002 and 2005, and on the other hand, a huge majority like in France in 
2002. 

PACS numbers: 02.50.Ey, 05.40.2a, 89.65.2s 



1 Introduction 

In the last years the study of opinion dynamics has attracted a growing number 
of works PElElElElEllZHHlinHinHnilI21 making it a major current trend 
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of sociophysics ITl] . Most models consider 2-state opinion agents combined 
with some local opinion update rule which implements the dynamics. They are 
found to lead to an opinion polarization of the whole population along one of 
the two competing opinions. A unifying frame was proposed to incorporate all 
these models JHl- Continuous opinion models yield similar tendency jlfil 117) . 

More recently, the concept of contrarian was introduced to account for some 
peculiar behavior of agents A contrarian is undistinguishable to others, 

i.e. its opinion evolves also by local rule updates. However, once it leaves the 
update group, it changes individually its opinion to the other one. The shift is 
independent of the opinion itself. A contrarian is not a permanent individual 
state. Each agent has a probability a to behave like a contrarian and (1 — a) to 
stick to its opinion. After each cycle of local updates, on average a proportion a 
of agents shifts spontaneously their opinion to the other one, while a proportion 
(1 — a) sticks to its current opinion. Given a fixed density of contrarians a, the 
associated opinion dynamics is then studied |18| . 

As intuitively expected, the existence of contrarians was shown to avoid total 
opinion polarization with the creation of stable attractors characterized by a 
stable coexistence between a large majority and a small minority. However, 
above some low density, they were found to produce an unexpected reversal of 
the dynamics, with the merger of the two attractors at the former separator, 
turning it to the unique stable attractor ^| Accordingly, for whatever 

initial conditions the dynamics leads to an exact global balance between the two 
competing opinions. It thus offers a possible explanation to recently observed 
hung election scenario |18| . A hung election being a two candidate run for which 
the result is very tied around fifty percent. Chaotic regime was also found in 
the description of investors in stock markets [231 • 

But in today campaigns, polls are regularly publicized making agents aware 
of which opinion is currently leading at the global level. It is therefore more 
natural to link the propensity to a contrarian behavior to the current level of 
global support for a given opinion. Accordingly, in this paper we relax above 
constraints of fixed independent contrarian density a in two successive steps. 
First, we study the effect on the dynamics making the density a proportional 
for each opinion to its current global support. Contrarians become propor- 
tionate contrarians. Second, we push the asymmetry further by restricting the 
contrarian behavior to only the current majority opinion, contrarians being one- 
sided. While proportionate contrarians are found to preserve the mean features 
of the fixed density contrarian dynamics, a more complex situation including a 
chaotic regime is discovered for one-sided contrarians. 

For a density a < a c — 1/9, one-sided contrarians produce a chaotic basin 
located in the fifty-fifty percent region. The associated Lyapunov exponent 
is calculated. However, there still exist majority-minority coexistence point 
attractors located on each side of the chaotic basin. Initial conditions determines 
which regime will dominate, either chaotic outcome around fifty percent, or 
point attractor with a well defined majority. For 1/9 < a < 0.301 only the 
chaotic basin survives. Further, in the range a > 0.301. the chaotic basin 
disappears and the majority starts to alternate between the two opinions with 
a staggered flow towards two majority-minority point attractors. 

On this basis the effect of additional social parameters are studied. A con- 
stant shift for contrarian activation is considered. Decoupling the local update 
time sequence from the contrarian behavior, and a delay in accounting for a 
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change of majority side are also introduced. Last, unsettled agents are included. 
They are found to only shrink the chaotic basin making the outcome very tied. 

Given above framework, the results may shed light to recent very unusual 
hung elections like in the 2000 US presidential vote and in Germany for the 
2005 elections. In both cases the outcomes were very tied with almost identical 
support. At the same time, the model can also provides an explanation for the 
2002 French presidential election where the winner obtained a huge majority 
around 80%. 

At this stage it is worth to stress that we are able to embody these contra- 
dictory voting outcomes within a single frame. But at the same time, we are 
not able to decide before hand which one will prevail. To overpass this difficulty 
would require the collaboration with social scientists to estimate the actual type 
of contrarians involved in a given election. 

The rest of the paper is organized as follows. In the next Section we review 
the original Galam model of contrarians. Throughout the paper, the size of 
update local groups is kept equal to three agents. The contrarians are then 
made proportionate in Section 3. Section 4 considers one-sided contrarians and 
the corresponding complex dynamics topology. Novel features are obtained. 
In Section 5 we study the effect of decoupling the local update time sequence 
from the contrarian behavior activation. A smoothing of the majority rule is 
also included. The effect of a status quo driven bias for contrarian activation 
is investigated in Section 6. Section 7 deals with the case of unsettled agents. 
Last section contains some discussion. 

2 The original Galam model: individual con- 
trarians 

We start recalling the Galam model of 2-state opinion dynamics model extended 
to the presence of contrarians ^S]- It considers a population of N agents where 
at a time t, NA(t) persons support one opinion A and {N — NA(t)} persons 
support another competing opinion B. In terms of global proportions among 
the whole population, it yields pt = for A and {1 — pt} for B. These 

values can be evaluated at any time using polls. 

From an initial value p t at time t, a dynamics is implemented in two steps. 
First, a neighborhood step where agents are distributed randomly among various 
size groups in which they update their respective individual opinion following 
the local initial majority 9 . In case of a tie in even size groups, either one 
opinion is adopted according to some probabilities 21 . The step is accounted 
by a discrete time increment of +1 leading to a new proportion pt+i of agents 
supporting opinion A. The second step is contrarian, each agent individually 
either shifts its respective opinion to the other one with a probability a, or 
preserves its current opinion with probability (1 — a). This second step yields 
an additional increment of time +1 and modifies pt+i to another value pt+2- 

Throughout this paper we explicit the calculations for the case of local groups 
with the unique value 3. Above rules thus yield respectively for the first step 
p t -> pt+i = Pmipt) with 

Pt+ i=P m (Pt)=P 3 t+3p 2 t (l- Pt ), (1) 
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and for the contrarian second step pt+i — * Pt+2 = Pc{pt+i) where 

Pt+2 = Pcipt+i) = (1 - a)p t+ i + o[l - pt+i}. (2) 

It should be stressed that before performing another cycle of opinion updates, 
agents are reshuffled [22] • Applying above 2-step cycle repeatedly n times re- 
sults in a proportion pt+2n of agents supporting A. All possible variations and 
extensions of the model can be included within a unifying frame |15| . 

At this stage we make a change of variable from p to d with p = d + i. It 
will appear to be more convenient for our investigation. A positive d makes A 
the majority opinion while a negative value grounds it as minority with a deficit 
|d| of support with respect to B. In terms of the new variable d, P m (pt) and 
P c (pt+\) become respectively 

dt+i = D m (d t ) = -2d? + |dt, (3) 

and 

d t+2 =D c {d t+ i) = (l-2a)dt+i, (4) 

which combine for one full cycle into the single Equation dt+2 — D c [D m {dt)], 
which we denote by 

d t +2 = D 2 (dt) = (1 - 2a){-2d 3 t + |dt}, (5) 

where index 2 of D% means one local rule followed by one contrarian step with a 
time interval of 2 for the dynamics. At this stage such a 2-step split could appear 
artificial but it will become instrumental latter on to extend the dynamics to 
cases where one cycle is built out of (k — 1) consecutive local updates followed 
by one contrarian effect. These cases will be denoted by d t +k — Dk{dt) with 
k being the appropriate time interval to study the associated properties of the 
dynamics. In the original Galam work k = 2. 

For this last case the main results obtained from Eqs. (1,2) or (3,4) are 
twofold. At low concentration a of contrarian behavior, total polarization is 
prevented with two mixed attractors at which a majority and a small minority 
coexist. The threshold for A victory is at p v = ^ or d v — which defines the 
separator of the dynamics. Furthermore, increasing a provokes a continuous 
phase transition at a c = | turning p v = ^ or d v = into the unique and stable 
attractor of the dynamics. The final state is a perfect equality of both opinions 
at the collective level with on going individual opinion shifts [T%| . 

3 Making contrarian behavior opinion current 
status dependent: proportionate contrarians 

While the original model considers a constant and fixed proportion a of con- 
trarians, it seems more realistic to make it depend on the current state of the 
system, in particular due to the existence of published polls. We thus suppose 
that if at some specific time t all agents are informed of the actual value of d t , 
they react accordingly as contrarians with respect to d t at time (t + 1) leading 
to dt+i- But it thus becomes natural to make the contrarian behavior opinion 
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dependent. Not to the opinion itself but to its current level of support in the 
population. 

To distinguish the associated contrarians from previous ones, we call them 
proportionate contrarians. Therefore, agents sharing opinion A react to p t , i.e., 
dt while agents sharing opinion B react to (1 — pi), i.e., — d t . We denote these 
rates respectively a(d) and b(d). From symmetry we have b(d) — a(— d) since 
opinions are time reversal. 

For the time being, we keep our 2-step cycle, which implies a regular and 
periodic publication of polls. Such a constraint will be relaxed in a latter Section. 
On this basis the proportionate contrarian density a becomes a function of time 
through the variable dt- If no agent shares the opinion A, no one will react 
against it as a proportionate contrarian yielding the constraint a(d = —1/2) = 0. 
On the other extreme at d = 1/2, a(d = 1/2) = ao where ao is the proportionate 
contrarian maximal value which satisfies < ao < 1. 

Keeping in mind that now the proportionate contrarian density is not the 
same at a given time among agents sharing respectively opinion A and B, Eq. 
(4) becomes 

df+2 = D 2 {d t ) = (1 - a t+ i - b t +i)d t+ i + &t+1 - Qt+1 , (6) 

where a t +i and bt+i respectively mean a(dt+i) and b{dt+\), and dt+i is given 
by Eq. (3). 

From Eq.(6) we can extract several properties of the associated dynamics. 
First we note that d = is a fixed point if and only if a(d = 0) = b(d = 0). 
It means no agent is contrarian at perfect equality of opinions. Second we can 
evaluate its stability by studying the value of the associated eigenvalue A with 
respect to one. When d = is a fixed point, it is an attractor when A < 1, 
which in turn implies the condition 

6a(d = 0) + a'(d = 0) - b'(d = 0) > 1, (7) 

where the prime means a derivative with respect to d. Otherwise when 6a (d — 
0) + a'(d = 0) — b'(d — 0) < 1 (A > 1), the fixed point d = is a separator. 
A separator implies the existence of two attractors located respectively on both 
side at d > and d < with a stable coexistence of a majority and a minority. 

Using above results we can review few specific functional forms for the a 
d-dependence. We start with the linear dependence a = agp = ao(l/2 + d). 
For ao < 1/5, d v = is a separator and the associated attractors are located 
at d = ±^==fi. When ao > 1/5, d = becomes the unique attractor 
of the dynamics. In other words, the former Galam result is recovered with 
proportionate contrarians stabilizing a perfect equality between both opinions 
once their density is larger than a critical value, here 1/5 ^B]- Results are shown 
in Fig. (1). 

Considering a power law form a — aop 1 = ao(l/2 + d) 7 we get the condition 
a o > 3+2-y t° make d = the unique attractor. It is a separator with two mixed 
phases attractors when ao < . Taking 7 = 2 yields a square dependence 
with | for the critical value of ao. More proportionate contrarians are needed 
( fl o > f = 0.286) with respect to the linear case (ao > \ — 0.20) to produce the 
perfect balance of opinions. When d = is a separator, the two attractors are 
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Figure 1: Application D 2 given by Eq. for the linear and symmetric case 
where a(d) = a^p = ao(l/2 + d) and b(d) — a(—d). On the left side ao = 0.1 < 
1/5 with two attractors at coexistence of a majority and a minority. On the 
right side ao = 0.3 > 1/5 with one attractor at a perfect balance among both 
opinions. 



located at d = ±|y 1 + ^-(1 + \fl + 80q). On the opposite, taking a square 
root dependence with 7 = 1/2, yields ~ 0.177 for the critical value of ao 
making less proportionate contrarians needed to get the opinion balance as the 
unique attractor of the dynamics. The former Galam model considers a = est. 
giving a critical value of | ~ 0.167 which is the lowest value we can get. 

At this stage we can conclude that proportionate contrarians do not modify 
qualitatively the former constant contrarian density dynamics opinion. 

4 Restricting contrarian behavior to the current 
majority: one-sided contrarians 

On the basis of above results, we go back to the case of a constant density 
of contrarians, but now restricting the activation of the contrarian behavior 
to only the current majority opinion. It thus yields for d > the conditions 
a(d > 0) = a — est (with < a < 1) and b(d > 0) = 0, while for d < the 
conditions are a(d < 0) = and b(d < 0) = a = est, where we have assumed 
symmetric conditions for both opinions. It is worth to note that this symmetry 
induces at d = the condition a(d = 0) = bid = 0) = | both preserving 
the same total density a of contrarians and recovering here the former original 
individual contrarian case. We call these contrarians one-sided contrarians. The 
associated rule update writes: 

d t+2 = D Q - B (dt + i) = (1 - a)d t +\ - | sign[d t+ i], (8) 

where sign[x] = 1 if x > 0, — 1 if x < and if x = 0. 

Using Eq. J3J for dt+\ and noting that sign[dt+i] — sign[d t ], Eq. © writes: 

d t+2 = D 2 {d t ) = (1 - o)(| d t - 2d\) - I sign[d t ]. (9) 
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At second order in d t it yields: 

D 2 (dt) — \ndt - ~ sign[ek], (10) 

where A m = |(1 — a) is the maximal slope of the application D 2 . 

We can now study the associated fixed points and their stability. We first 
note that the point d = is a singular fixed point. Then, the application being 
symmetric we restrict the study to d > 0. According to Eq. 10, fixed points 
exist ifa<a c = i~0.11 and their values are: 

When they exist, aV_ is unstable and aV+ stable (see Fig. 2). From Eq. (fTTTf) 
it gives: 

= w^rr (12) 

with the eigenvalue A ~ A m > 1 (A m > 1 for a < a c ). 

Eq. I|10(l exhibits clearly condition for expansion when A TO > 1, and fold- 
ing up by the discontinuity at the origin d = 0. Thus, for being chaotic, the 
application D 2 has to satisfy two conditions; first, the expansion condition and 
second, the perpetuity of the chaotic basin [23 E]. If the application D 2 is 
chaotic, the interval of successive iterated points after transients is called f2;. 
Here fli =] - a/2; a/2[ (see Fig. 2). 

1 st condition : expansion. 

Inside the interval tti the application D 2 possesses a slope A, such as A m > 
A > D' 2 (a/2), where D' 2 denotes the derivative of D 2 relating to d. This implies 
that necessary A m > 1, i.e. a < 1/3. Furthermore, D' 2 (a/2) = 1 for a ~ 0.27. 
According to this unique condition, D 2 loses its chaotic nature for a including 
to the two last values. 

2 nd condition : perpetuity of the chaotic basin, i.e. D 2 (p,i) C Hi- 
If there are not fixed points, this condition is automatically satisfied. If fixed 
points exist, i.e. for a < a c , they must be outside the interval ] — a/2; a/2[, i.e. 
dp- > a/2. According to Eq. (|12|) it implies A m < 2. This condition is always 
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Figure 3: Successive iterated points by the application Di given by Eq. 10. 
Notation here: d(n + 1) = D% \d{n)\ . Left side: a = 0.05 with the initial value 
d(0) = 0.05. Right side: a = 0.2 with the initial value d(0) =0.1. 
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Figure 4: Error growth and Lyapunov exponent of the application D2 given by 
Eq. © in semi-log plot. Notation here: d(n + 1) = D2[d(n)l . The initial error 
is <5(0) = 10~ 10 , taken after transients. Left side: a = 0.05. From the graph the 
Lyapunov exponent is Xi yap — 0.353 while ln(A m ) ~ 0.354. Right side: a = 0.2. 
From the graph the Lyapunov exponent is Xi ya p — 0.166 while ln(A m ) ~ 0.182. 



satisfied here for group size 3 while it is not the case for larger update groups 
from size 5. Nevertheless, to observe a chaotic behavior initial condition must 
satisfy \d[t = 0)| < dp-- Otherwise, the successive iterated points go to a stable 
fixed point at ±dp+. 

A chaotic behavior is numerically observed until a ~ 0.301. Illustrations are 
shown in Fig. (3) for a = 0.2 and a = 0.05 (with an appropriate initial condition 
in the latter case). 

Fig. (4) shows sensitivity to initial conditions and permits the evaluation of 
the Lyapunov exponent Xi yap - The initial difference 8(0) grows exponentially 
5(n) ~ 6(0) e nXlyap , until saturation at the typical size of the interval il/. The 
Lyapunov exponent is positive with \i yap — ln(A m ) as seen from Eq. (|10Jl . 

The notation adopted in these figures is to consider the n th iteration of the 
application D2 since the beginning of the electoral campaign as the equivalent 
time n. Thus dt +2n = d(n), where to is the time at the beginning of the 
electoral campaign. In other words, d(n + 1) = D% [d(n) \ . 

At this stage we can make the following comments: 



8 



0.02 0.04 0.06 0.08 0.1 

a 



Figure 5: Separators ±dp- of chaotic basin in dashed line and point attractors 
±aV+ in plain line, given by Eq (|11J) as functions of density a. a < a c = i. 

- The discontinuity of the application Di at d = is not the unique 
origin of its chaotical nature. Indeed, we can transform Di to be continuous 
and derivable, e.g. via a — > a(d) — a m (1 — e~~). If a m < 0.301 and a <C 1 then 
this application exhibits a chaotic behavior. Let us recall that sudden variation 
is the source of the chaotic nature of the application because it permits the 
folding up of an expansive application (if A m > 1). 

- When a > 0.311 the application has stable fixed points of doubling 
period. They are obtained from D2(d) = —d. 

- Increasing a, the first separation for d > inside the interval Qi in 
two intervals occurs at a ~ 0.056. This result can be retrieved considering 
the doubling iterated application = D2 o D 2 , i.e. D 2 2 \d) ~ /^(-^(d))- 
Indeed, this occurs when lim D 2 2 \d) > d (2) *, i.e. D 2 (a/2) > d {2) * , where 

d^* is a fixed point of doubling period where d^ 1 * is obtained via D2(d^*) — 
-d^*. From Eq. UHJ) we retrieve a ~ 1 - 2^1 ~ 0.057. The extension 
of successive iterated points after transients, Qi, is then, for positive values: 
]0;C 2 (a/2)[U]|^ 2) (a/2)| ;a /2[. 

- Starting for instance from opinion A being initially the majority, i.e 
d > 0, we evaluate d c h, such as D2(d c h)=0. Thus, if d < d c h then D2{d) < 
and the majority side will change to the opposite side; and reciprocally if d > d c h 
then D2(d) > and the majority side will keep the same side (see Fig. 2). From 
Eq. dJ: 

dch =s (13) 

Since successive iterated points are contained into the interval Qi =] — a/2; a/2[ 
and I?2(i/2) < d c h with X m < 1.5, we deduce that d > cannot remain 
positive more than twice before turning negative. Note that at second or- 
der on d, 1)2(0/2) > d c h yields a second order equation whose solution is 
Am > ^ 1.62, the golden number. 

To summarize, for a < a c — ^ this model provides the coexistence of three 
radically different basins. A chaotic one located around d = delimited by 
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separators ±dp_ and two others at the extremes with each one a point attractor 
±dp+- Fig. (5) shows positions of separators and points attractors as functions 
of a. For 1/9 < a < 0.301 only the chaotic basin survives. For a > 0.301 the 
opinion forming staggers towards two stable fixed points of doubling period. 
Most interesting effects are obtained for a density of one-sided contrarians such 
that a < a c = |. Indeed, in this case the initial value e?(0) at the beginning of 
the electoral campaign determines which one basin is selected. If |rf(0)| < dp _ 
the intention vote dynamics is chaotic but with a result around d = such that 
d e] — a/2;a/2[. If \d(0)\ > dp-, the issue is certain but with a result at the 
extremes. This could account for contradictory electoral outcomes whose results 
are either around p = 50% or with a huge majority at p ~ 80% like the 2002 
French presidential election. 

5 Varying the time scales of collective informa- 
tion and local updates 

After having considered 2-step processes, we now study the opinion dynamics 
driven by a k-step process where the individual activation of the one-sided con- 
trarian step occurs after k — 1 repeated steps of local majority rule update. It 
accounts for the fact that polls are not made public every single day during a 
campaign while on the contrary people keep on discussing all the time. 

Now d t +k = D _ s (d t +k-i) with d t+ k-i = D m l ^ 1 \d t ) where Dm ^ = D m o 
D m o ■ ■ ■ o D m , i.e. k — 1 iterations of D m . In previous sections, k was equal to 
2. Accordingly one-sided contrarians consider the collective information with a 
delay acting at time t + k — 1 while considering information at time t. However 
they could as well consider the collective information without delay, just after 
the last update inside groups at time t + k — 1. Indeed, they act according to 
sign[d) (see Section 4) and sign[d t ] = sign[dt+i] — ■ ■ ■ — sign[dt+k-i)- 

Moreover dt+k-i increases very quickly from the origin d = 0, making a 
contrarian behavior without much effect. For instance k = 10 makes the slope 
at the origin at (|) 9 ~ 38. So, we have to slow down the dynamics driven 
by the update rule inside groups in order to study the varying time scale of 
collective information. It is done quite naturally assuming that not every agent 
eventually changes its opinion to follow the local majority within each cycle of 
local updates turning Eq. (JIJ to: 

Pt+i = P m , w {pt) = w[Zp 2 t - 2p 3 t ] + (1 - w)pt, (14) 

where w denotes the propensity of an agent to be convinced by majority rule 
with < w < 1. Using the d variable gives: 

d t+1 = D mtW (dt) = (1 + w/2)d t - 2w d 3 t . (15) 

Now k — 1 iterations of D m w yields a slower slope at the origin, e.g. for 
w = 0.1 and k = 10 it is (1 + w/2) 9 ~ 1.5. 

Let a the density of contrarians and w the propensity of an agent to be 
convinced by majority rule. The new intention vote dynamics of the k-step 
process, generated by the one-sided contrarian step occurred after k— 1 repeated 
local majority rules writes as 

d t+k = D k (d t ) = D - s [D^(d t )] , (16) 
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Figure 6: Successive iterated points by the application Dk given by Eq. (|TB|l . 
Notation here: d(n + 1) = Dk[d(n)\. Left side: a = w = 0.1, fc = 3 and an 
initial value d(0) = 0.1. The application is not chaotic. A m ~ 0.99 < 1. Right 
side: a = w = 0.1, fc = 18 and an initial value d(0) — —0.04. The application is 
not chaotic by escaping from the previous chaotic basin. A m ~ 2.06 > 2. 



where Z? _ s and D m _ w are respectively given by Eqs. (0 E3>- This yields at 
second order on dt the Eq. (|1C)|> . but now with the slope A m = (1 — a)(l + ^) k ~ 1 . 
To exhibit a chaotic behavior the condition for expansion A m > 1 gives now: 

fc-i> -y 1 "! (17) 

ln(l + w/2) v ; 

If a, w 1 and are of the same order, then fc — 1 > e.g. for a = w <C 1, 
fc > 3. Numerically this can be satisfied until a = w < 0.4 (see Fig. 6). 

With respect of the perpetuity of the chaotic basin, the unstable fixed points 
(Eq. H12fl ) are idp- ~ ± 2 (\ °-i) ' ^ they exist. It implies A m > 1 and A m < 

1 + a to have di?- < 1/2, so, fc — 1 > g~n+W2) • ^ °' w ^ •"■ anc ^ are on ^ ne same 
order, then fc — 1 > — , e.g. for a = w <C 1, fc > 5. Nevertheless, the second 
order approximation is not proper as soon as the inequality dp- < 1/2 doesn't 
satisfy \d\ <C 1. 

To check if successive iterated points cannot escape the chaotic basin, i.e. if 
dp- > a/2, from Eq. I|12(l we need to have X m < 2, i.e. 

*-i< 'T»- h V ) - (18 » 

ln(l + w/2) 

If a, b <C 1 and are on the same order, then fc — 1 < 21n( ^ +2a , e.g. for 
a = w = 0.1, fc < 16.9 while numerically, fc < 18 (see Fig. 6). As in Section 
4, initial value has to satisfy \d(t — 0)| < dp-- However, contrarily to the 
precedent section, here successive iterated points can escape the previous chaotic 
basin. 

In addition, for fc > 2, this k-process with fc — 1 repeated steps of local 
majority rule updates increase naturally the majority side persistency before 
changing compared to the previous model. 
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6 Status quo driven bias 



Up to now both opinions were perfectly symmetric. However while dealing with 
political opinion dynamics in view of an election a difference should be made 
between the opinion supporting the current political party in power and the one 
supporting the challenging party. The contrarian behavior should be a little 
more active against the former winner opinion thus creating a bias in favor of 
the challenging opinion 9 . 

Assuming B is the former opinion winner and denoting s the bias in favor 
of the challenging opinion, Eqs. I10|) should be rewritten as: 

d t+2 = D 2 (d t ) = (1 - o)(| d t - 2d?) - ~ sign[d t - a], (19) 

and Di(dt) — \ m dt — § sign[dt — s] where A m = f (1 — a). 

In the case a < a c = 1/9 for which fixed points exist, Fig. (7) (a = 0.05 and 
s = 0.02) shows that for \s\ < dp-, the bias doesn't modify the position of the 
fixed points. But now, depending on the ration s to a, the successive iterated 
points could escape from the previous chaotic basin, and thus reach a point 
attractor. The values s, at a fixed density a, for which this new phenomenon 
occurs, is obtained when Sup( lim Z^^t) ; lim |£>2(dt)|) > dp-. 

dt — >s~ dt — >s+ 

Thus, at first order on s, (a/2 + A m |s|) > dp- and at second order on dt, 
Eq. lO yields: 

\s\>s c ^a 2 ~ Xm (20) 

where here 1 < A m < §. For instance a = 0.05 gives s c ~ 2.37% for the exact 
value s c ~ 2.44%. 

Evaluating the values s, at a fixed density a, for which successive iterated 
points dt have the same sign we find \s\ > d c h, where d c h is defined as D2(d c h) = 
for the application without bias (see Fig. 7). From Eq. (|13fl : 

\s\>a^-. (21) 

For instance a = 0.05 yields \s\ > 1.75% for the exact value \s\ > 1.76%. It is 
worth noting that for s sufficiently large the application is no more chaotic. For 
instance, with a = s = 0.2, the application has a periodic attractor of period 
13. 

Including a bias in favor of one opinion provides two main effects. First 
one, after transients, a majority side (in the favored opinion) before changing 
more persistent than without bias. Second one, for given density a of one-sided 
contrarians, it is possible to escape the chaotic basin without bias to reach the 
point attractor of the favored opinion. 



7 Unsettled people and contrarian attraction 

We now go back to a population where agents sharing either one opinion evolve 
by local majority rule updates only without contrarians 9 . But we also consider 
another population of agents who do not take part in the public debate. They 
are unsettled and hold no opinion. However they are gradually driven in the 
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Figure 7: Application D2 given by Eq. (|19fl with a status quo driven bias for 
a = 0.05 and s = 0.02 in plain line. The dashed line is the application without 
bias. At the right side, a zoom of this application around d = 0. 



public debate on a contrarian basis. At a constant rate u with < u < 1 they 
move to the opinion holding population starting with an opinion opposite to the 
current majority. Once they adopt an opinion they become identical to other 
sharing opinion agents, i.e., they evolve by local majority rules. 

At time t the number of persons sharing opinion A, opinion B and no opinion 
are denoted respectively N A (t), N B (t) and Nu(t) with N (t) = N A (t) + N B (t) 
and No(t) + Njj(t) = N where N is the total number of agents of both popu- 
lations. Associated probabilities are: 

N A (t) Nb{^ 
Pt = — , , and » T . . = 1 — p t ■ (22) 
N {t) N (t) Pt v ' 

We still have a two-step process. The first one is unchanged with pt+i = 
Pmipt) while the second one is produced by the contrarian unsettled agents 
partial joining the debate. Note that now the application p t — > pt+2 is no longer 
stationary. To account for the shrinking dynamics of unsettled agents we note 
n the time which corresponds to the n th iteration i.e. dt +2n = din) where to is 
the time at the beginning of the campaign. Thus, writing d(n+l) = D2. n > 
it gives: 

d(n + 1) = D 2 , n [d(n)] = * d( " } ~ ^ a } H ) Si gn[d(n)}, (23) 

l + a(n) 2(l + a(n)j 

with 

, \ N u(n) (l-u) n 
a(n) = u = u— ^ — !—— , (24) 

V ; N (n) R- (1-u)"' V ' 

where R = > 1. For u < 1, on the first order on Eqs. EJJ 

are unchanged with now instead of a an effective time dependent contrarian 
density: a — > a(n). 

In the limit u <C 1 from Eq. Ijl2|l we can write, d nt F- — a(n) with 0; jn 
defined for < |d| < With lim a(n) = the interval „ shrinks to 0. 

n — >oo 

However this model prohibits n — > 00 since the number of unsettled persons 
mobilized at each iteration is an integer and not a real one, although the n th 
iteration, the number uNu(0)(l — u) n is assumed to be greater than 1. 
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Figure 8: Applications Z?2,ri given by Eq. I|23() with u = 0.02 and R = 2. n = 
in plain line and n = 1 in dashed line. 

In the case (R - 1) > it, a(rt + 1) ~ (l - u -^-) a(n) from Eq. pijl. 
Accordingly points escape the basin bordered by ±d nt F- from time nton+1, 
if \d(n + 1)| > d n +i,F- — a(n + 1) implying \d(n)\ > (1 - u 2 (R R i) ) a ( n ) ( sec 
Fig. 8). Thus, if d belongs to the interval fij iW at time n, successive iterated 
points will be contained into the successive intervals f2z, ra . (see Fig. 9). 

Even if successive iterated points do not escape the basin bordered by 
±d ni F-, the dynamics is no longer chaotic. Indeed the point d = is now 
asymptotically stable. This is due to the non-stationary dynamics effects. Nev- 
ertheless, the shrinking dynamics doesn't affect the sign of iterated d. So, the 
expected winner at the issue of an electoral campaign remains unpredictable. 

Including unsettled people driven gradually in the public debate on a contrar- 
ian basis provides a shrinking of the chaotic basin of the one-sided contrarians 
model. Thus, the results may shed light to recent very unusual elections like the 
hung 2000 USA presidential and German 2005 elections. Indeed, the winners of 
these long electoral campaigns were very unpredictable and the outcomes very 
tied. 

8 Conclusion 

We have presented a simple model giving a deterministic opinion dynamics (see 
Eq. EJ, which can be chaotic. Furthermore, chaotic basin can coexist with two 
point attractors at the extremes. This model contains two main effects. First 
effect, amplification one given by the local majority rule inside groups. Second 
effect, retroaction one given by the action of the one-sided contrarians. The one- 
sided contrarians act by comparison and opposition to a collective information, 
the majority. Their action introduces to dynamics a discontinuity. 

Afterwards, rooted in this simple model, some others features of electoral 
campaigns are added, like the fact that polls are not made public every day, 
or a bias in favor of one opinion, or the influence of unsettled people gradually 
driven in the public debate on a contrarian basis. 

To sum up our various results we come back to the one-sided contrarian 
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Figure 9: Successive iterated points with D 2y n given by Eq. l|2*3)l where u = 
0.02 and R = 2. Borders ±^S^ are included in thin points. Notation here: 
d(n + 1) = D 2 , n [d(n)] . 




model deterministic equation used throughout this paper: 

/Or) = (1 - a)(~x - 2x 3 ) - | sign[x], (25) 

where x = d G [—1/2; 1/2] and a is the density of one-sided contrarians, < 
a < 1. At the second order on x approximation, with A m = |(1 — a), fix) ~ 
A m x — ^ si<?n[:z;]. A m is the expansion effect. The term — | sign[x] is related 
to the folding up, i.e. the retroaction effect. The term — x 3 can be seen as a 
non-linear saturation effect and gives the point attractors. It is symmetric with 
respect to x = 0. 

Afterwards, accounting for some others social parameters enrich the model 
and modify a little bit its applications. By decoupling the local update time 
sequence from the one-sided contrarian behavior activation, the two effects, 
amplification and folding up effects, can act separately. Next, a bias in favor of 
one opinion introduced as a simple parameter s fixes the discontinuity position at 
x — s, i.e., sign[x] — > sign[x — s]. The opinion dynamics is no more symmetric. 
Last, unsettled people driven gradually to public debate on a contrarian basis 
roughly plays with the folding up effect uniquely; A m ~ | and a — > a(n) with 
a(n) — > when n — + oo. The opinion forming dynamics is no more stationary. 

At this stage it is worth to stress that contrary to what could be expected, 
our treatment using probabilities, local updates and reshuffling between updates 
does not define a mean field like frame. This result was demonstrate recently 
using Monte Carlo simulations of a nearest neighbor ferromagnetic Ising system 
on a square lattice |22| . Indeed it creates a new class of universality in addition 
to both 2-d Ising and mean field ones. For our current model, in the case of 
group size four with no contrarians, using a cellular automata was shown to 
recover our analytical result |12j . 

To conclude it is worth to notice that in terms of real life situations, due to 
the finite number of iterations imposed by the fact that any public campaign 
is finite in time, the intention vote dynamics will not exhibit a chaotic behav- 
ior although in principle it could. In addition, the non zero fuzziness of poll 
measurements of the initial intention vote distribution, result automatically in 
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Figure 10: Combination of a k-step process and unsettled agents with k = 10, 
w = 0.05, u = 0.05 and R = 1.5. Initial values are 0.1 in plain line and 
(0.1 - 0.001) in dashed line. 




a growing error making difficult to predict the sign of successive iterations (see 
Fig. 10). 

Although this simple model does not pretend to account for an exhaustive 
explanation of opinion forming during electoral campaigns it exhibits some fea- 
tures which could shed new light on recent surprising voting outcomes, like for 
instance, on the one hand, an unpredictable issue with a very tied outcome like 
the German 2005 and the USA 2000 elections, and on the other hand, a well 
predicted outcome with a huge majority like in the 2002 French presidential 
elections with a majority around 80%. 

However it is worth to stress that at this stage if we are able to embody 
above contradictory voting outcomes within a single frame, we are not in a 
position to select the prevailing one. To overpass this difficulty would require 
the collaboration with social scientists to estimate the actual type of contrarians 
involved in a given election. It is open for future work. 
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